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Abstract 

A new method is developed to calculate multiloop Feynman diagrams in 
non-relativistic field theories. A consistent scheme for regularizing and renor- 
malizing loop integrals is established and shown to reproduce the results of 
dimensional regularization (DR) and modified minimal subtraction {MS) or 
DR and power divergence subtraction (PDS) up to next-to-leading order 
(NLO). However, significantly less effort is required to evaluate the integrals 
and the methods are easily generalized to higher order graphs. Thus, even the 
most complicated multiloop graph can be expressed in terms of analytic func- 
tions. These techniques are then used to calculate the quadrupole moment 
of the deuteron to three loops in the Kaplan, Savage, and Wise effective field 
theory. A new unfixed direct S^D wave counterterm occurs at this order 
and its value is determined. 
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I. INTRODUCTION 

Recent progress in effective field theory Q has made possible the cal- 
culation of many properties of nuclear systems in a perturbative expansion 
|2| ||]. The progress is largely based on a consistent power counting developed 
by Kaplan, Savage, and Wise (KSW) which accounts for the large scattering 
lengths in the two nucleon sector by utilizing dimensional regularization and a 
unique subtraction scheme called power divergence subtraction (PDS). How- 
ever, higher order graphs very quickly become intractable with the standard 
techniques of dimensional regularization and Feynman parameters. This is 
because combining three dimensional propagators in the usual way leads to 
square roots of polynomials of the Feynman parameters and one quickly be- 
comes entangled in a web of hopeless integrals. While these may be computed 
numerically for finite loop graphs, many of the more complicated divergent 
loop integrals remain elusive. Thus, one would like to find a new method for 
evaluating graphs which avoids the shortcomings of these techniques and also 
allows for consistent regularization and renormalization schemes. Further- 
more, it would be beneficial to maintain the desirable KSW power counting. 
Other subtraction schemes |7L|j have recently been developed that obey the 
KSW scaling, but none of these simplify the calculation of higher order con- 
tributions. 

In this spirit, a consistent regularization and subtraction scheme is de- 
veloped for non-relativistic field theories. In particular, the KSW effective 
field theory is used throughout to develop the methods. The same recipe 
is applicable to all diagrams. Loop integrals are transformed to coordinate 
space and divergences are transformed into divergences in a parameter with 
units of energy, where they are consistently treated. The results are seen to 
be totally consistent with dimensional regularization and PDS (or MS) in 
the known cases. In addition, the coordinate space-parameter subtraction 
method, which we choose to call CPS, is easily extended well beyond the 
realm of the standard techniques in calculating higher order loop graphs. In 
section II the general methods are developed for bound state graphs. We show 
how calculate 2-pion exchange bound state graphs with an incident photon 
in the zero momentum limit, which is necessary for the NNLO calculation of 
the deuteron quadrupole moment. In section III, the results are generalized 
to finite photon momentum transfer graphs and above threshold scattering 
graphs. These examples serve to illustrate how graphs of any order may be 
obtained in terms of analytic functions. Section IV reviews the KSW effective 
field theory and the calculation of electromagnetic properties. Section V de- 
tails the calculation of the deuteron quadrupole moment and the results are 
shown. 

II. BOUND STATE GRAPHS 

In Ref.[5], the deuteron electromagnetic form factors are calculated to 
NLO. Calculation of the full form factors at NNLO requires several graphs, 
the most difficult of which is the 2 pion exchange with a finite momentum 
photon attached on one side (Fig.l). First, we show how to evaluate the 
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FIGURES 




Fig.l Two potential pion exchange graph that gives a contribution to the deuteron quadrupole 
moment at NNLO. The photon corresponds to A . 

graph at zero momentum transfer. This will serve to introduce the CPS 
method, since this complicated graph contains almost all possible classes of 
loop integrals that might occur in bound state problems. 

In the EFT of Kaplan, Savage, and Wise we find an amplitude for Fig.l 

of 

_ -9g A M%e t dpdldk(2l i k j l-k - kljk 2 - hkjl 2 + \l 2 k 2 5ij) 

~ 2/4 J ((p - q/2) 2 + 7 2 )(p 2 + a 2 )((p + I) 2 + a 2 )((p + l + k) 2 + a 2 )(Z 2 + m 2 )(fc 2 + m 2 ) 

(1) 

where dp = ~A^s , Q 1S the photon 3-momenta, p =p , I = I , k =k , 7 = 

^—EMn, = 138MeV is the pion mass, a = \fj 2 ~+~^-, and i,j are the 
initial and final deuteron spin indices, respectively. We wish to extract the 
quadrupole moment part at zero momentum transfer. After straightforward 
but tedious convolutions and projections, the numerator of eq.(l) can be 
written in terms of p 2 , I 2 , k 2 , (p — I) 2 , and (p + k) 2 as 54 independent terms 
such as (p — l) 4 (p + k) 2 l 2 k 2 /p 4 , k 2 l 4 , and other products of squares of the 
momentum variables. Note the appearance of l/p A terms; these will introduce 
artificial infrared divergences which must cancel in the end. After all these 
steps, one is left with integrals over products of 1, 3, 4, and 5 propagators. 
Note that the (p — l) 2 propagator is raised to the fourth power from the Taylor 
expansion. 

Let I(n), J(n), K(n), and L(n) be 1, 3, 4, and 5 propagator integrals with 
the (p — I) 2 propagator raised to the nth power. We now explicitly show how 
to calculate these integrals in CPS. 

1-PROPAGATOR INTEGRALS 

By definition 7 7 (n)= / dp -^^ryi ■ Transforming to coordinate space we 
have 
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/ 7 (1) = J dp J d 3 xG 7 (x)e-^'~ i 



where G 1 (x) = is the Fourier transform of the propagator. Trivially 

I 7 (l) = G 7 (0) = 4^(1/^ — 7)U=o- Labelling the divergence /i we have 

Alternatively we could have differentiated with respect to 7 and integrated 
back, yielding 7 7 (1) = — f^ d / y2 , yfd 3 x(G- / (x)) 2 , which becomes / 7 (1) = 
~~ Joo^T^F = 7 ^ XT 7 , Labelling the linear parameter divergence as \i we have 
eq.(3). In fact, eq.(3) is exactly what would be obtained using dimensional 
regularization and PDS. This is a general result : poles in 3 dimensions 
using dimensional regularization and PDS are linear parameter divergences 
in CPS. Other I 7 (n) for n > 1 are obtained easily by differentiating eq.(3) 
with respect to 7. We also have integrals of the form IO^fdp-ij, for which 
we let 7^0 in the above. Hence, 

An 

The standard DR dictum that scaleless integrals are zero seems to indicate 
that CPS is not equivalent to dimensional regularization and PDS in eq.(4). 
However, scaleless integrals should only be zero when 4 dimensional logarith- 
mic divergences are considered, such as in MS, since then there is only one 
scale (the renormalization scale) with which to make a dimensionless argu- 
ment for the logarithm. But for power law divergences, as in PDS, all we 
need is one scale, the renormalization scale. Thus, the natural extension of 
dimensional regularization to PDS is equivalent to eq.(4). Finally, we consider 
IZ=fdp-^. After Fourier transforming and performing the integrals, we are 
left with IZ = ^L j^dx. Labelling the linear coordinate divergence as v we 
have 

IZ = ^. 

The same result is obtained by letting 7^0 in 1(2) = and then ^\ 1= q^2u. 
For our purposes it is not really important how one labels these artificial 
infrared divergences since a conspiracy amongst the 54 terms must cancel 
them in the end. 

3-PROPAGATOR INTEGRALS 
By definition 

J(n)=Jdpdl — 2 - -^j- _ l)2 - ^ 2)n(/2 - . 

In coordinate space we have 

J(l) = Jdpdld 3 xd 3 yd 3 zG 1 (x)G 1 (y)G m ^{z)e- i ^^ + ^- l >y +l - 7 l 



This quickly reduces to J(l) = Jd 3 xd 3 yd 3 zG 1 (x)G 1 (y)G m7V (z)5(x + ^y)5(y 
~ z ) = 16^ f(?dx^-^-, where b = 27 + m,,-. Differentiating with respect to b, 
we have J(l) = j^rf^dbl = ^[logft - log6|oo]- Hence, 

Tf1\ ~ l 1 2 7 + "V 
J(1) = 16^ l0g ^^- 



Again, this is the same result obtained using DR and PDS or MS, but 
those require substantially more effort. In general, four dimensional logarith- 
mic divergences of dimensional regularization become logarithmic parameter 
divergences in CPS. Using the same methods we can derive 

JO(1) 55 l dpdl p*{( P -l?+l*W+ml) 



167T 2 fjL 



and 



JZ(1) fdpdl p4({p _ l)2 + ^ ){p + m2) 

_ 2u(m n + 7) - 1 
~ 32^(7 + ^^)2 ' 

where v is defined after eq.(5). Of course we get J(n), JO(n), and JZ(n) for 
n > 1 by differentiation with respect to 7. 

4-PR0PAGAT0R INTEGRALS 

There are several different classes of 4-propagator integrals. The first is 
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K{n)= / dpdldk 



' ((P - 2 + 7 2 ) n ((P + fc) 2 + 7 2 )^ 2 + m 2 )(k 2 + ml) ' 

which contains both 3 and 4 dimensional poles in dimensional regulariza- 
tion. In configuration space we obtain K(l) = jd 3 x{G 1 {x)) 2 {G mw {x)) 2 = 
^jj^dx^-^-, where a = 2(m n + 7). By partial integration K(l) = 

64^[ — ~F~lo° ~~ a /o°^ x ^~"]- Hence the divergences are manifested as the 
linear and logarithmic poles seen before in eq.(3) and eq.(8). The result is 

K(D = e^[2(m» + 7) Oog 2( ^ +7) - 1) +/*]■ 
Next we consider integrals of the form 

KP(n)= J dpdldk 



((P " 2 + 7 2 ) n ((P + ^) 2 + 7 2 )(/ 2 + m 2 n )(k 2 + mj) ' 



This is greatly complicated by the appearance of p 2 in the numerator which 
cannot be decoupled with momenta in the denominator. But one can evaluate 
the integral simply by using p 2 = Jd 3 x5(x)(—Vl)e~ tp ' x . It follows that 

KP(1) = -Jd 3 xd 3 yd 3 zG y (y)G 1 (z)G m7T (y)G m Jz)5(x)V 2 J(x + y + z ) 



where the derivative operator acts only on the delta function. Now we use 
Jd 3 xf(x)V 2 5(x) = f"\ x=0 , where the primes denote operation by V 2 . Thus, 

KP(1) = -J d 3 yd 3 zG^y)G^z)G m ^y)G m ^z)8"{y + 7) 
= -|d 3 yG 7 (y)G^(y)[G 7 (y)G m Jy)] // 
From the definition of G 1 {x) we find 
V,, G 7 (x) = 

V 2 x G 1 {x)= 1 2 G 1 {x)-5{x). 

Straightforwardly, 

KP{1) = G 7 (y)(G m „(y)) 2 \ y=0 + (G 7 (y)) 2 G mn (y)\ y=0 - ( 7 2 + m 2 ) J d 3 y(G y (y)) 2 (G m 
= I 2 {l)I m „ (1) + J 7 (1K1 (1) - (7 2 + ml )K(1). 

We must also deal with KPP(n), which is the same as KP(n) except the 
numerator is p A instead of p 2 . This is obtained in similar fashion and the only 
new result needed is V x G 1 {x) = 7 4 G 7 (x) — r y 2 5{x) — 5"(x). We obtain 

KPP{\) = ( 7 2 + ml) 2 K{\) - 2(2 7 2 + m 2 )/ 2 (l)/ m ^ (1) - 2( 7 2 + 2m 2 )/ 7 (l)/^ (1). 



5-PROPAGATOR INTEGRALS 

Finally, we come to the master integrals for 2 pion exchange graphs. By 
definition 
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L(n)= j dpdldk 



(p 2 + a 2 )((p - I) 2 + 7 2 )"((;p + k) 2 + 7 2 )(/ 2 + m 2 w )(k 2 + m 2 ) ' 



where a will be set to 7 in the end, but we label it differently for reasons 
that will be clear. After the delta function integrations in coordinate space 
we arrive at 



L(1) = Wfj d ' ixd6 y 



e -(~f+m n )(x+y)-a\ x-y\ 



(4vr) 5 7 x 2 y 2| x -y \ 

1 roo roc e -(n+m w ){x+y) 

"2a(47r)3io Jo V xy [& 1 
= 2a(4vr) 3/L ' 

where we have performed the angular integrations in the last step and defined 
II in the obvious manner. Now use the symmetry between x and y to write 



I L = 2 (°°— f X ^y e -(7+">-)(x+y){ e -a 
J x J y 



(x-y) _ e -o(x+j/)i 



- it - iE- 
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One must be careful here. We used the property that the integrand is sym- 
metric in x^y to get rid of the troublesome absolute value. But this must 
be done for both pieces, since L(l) is finite but each individual piece l£ and 
if is divergent. The wrong answer would be obtained by exploiting the sym- 
metry in l£, but not if. This is a consequence of the following simple rule : 
only finite integrals may be exploited for their symmetry properties as above, 
since otherwise information is lost about the relationship between the artificial 
divergences. In splitting a finite integral into several divergent pieces, every 
piece must be dealt with in the same manner. Now if = 2/~^ jxd^ e -bx-my ^ 
where b = m^+j+a and rh = m^+j — a. We only need evaluate l£ explicitly 
since if = l^(rh—^b). Differentiating with respect to b and rh we have 

rrh rOD Ary rx 

1* - -0 I A™ I —I ^,^bx-rhy 



/m rOD rx 

dm / — / dye' 
00 J x J 

/rh r Jrh rb roc 
™ db dxe- bx (e-™ x -1) 
00 m J 00 J 



= -2 / — / db(- 



00 m J 00 rh + b b 

„ f m dm m + b , b . 
= -2/ —(log log-), 

where we have used the CPS prescription for the artificial logarithmic diver- 
gences. Next we employ 



/ 



log (x -\- yj x 
dx = log x logy - Polylog 2 ( ), 

x y 



whence we deduce 

It = 2Polylog 2 (-^) 

The poly logarithms are generalized logarithms and may be defined by 
Polylog„(x) = Ya^i fs"- It can easily be shown that for large x, 

Polylog 2 (-xH - - T- Thus . 

,f = 2[Po ly log 2 (^) + ^ + £] 
and from this we obtain if = 2[ < -^- + fj], so that I L = 2[Polylog 2 (^) + 

2 

Y2]. Hence, we finally have 

L 1 = 7ji3- Pol y lQ g2 1 — 1 + 

1 r — m w 7r 2 



■[Polylog 2 ( 



(4tt) 3 7 L j fc52V m 7r + 27 / 12 J 



L(2) may be obtained by differentiating this with respect to 7. One might 
hope to get L(n) for n > 2 similarly, but because the x<^y symmetry was 



required, we can not label the 2 7's in eq.(19) differently and then differentiate. 
However, we can obtain L(3) and L(4) by differentiating eq.(26) with respect 
to 7 and using the following integral which is easily calculated : 



L22 = J dpdldk 



1 



(p 2 + a 2 )((p - I) 2 + 7 2 ) 2 ((p + k) 2 + 7 2 ) 2 (P + m 2 )(k 2 + m 2 ) 
1 



87 2 (47r) 3 (m 7r + a + j) 2 (m n + 7) ' 
By explicit calculation or by taking a— >0 in eq.(26) we arrive at 

1 



LO(l) = / rfpdfdfc 



P 2 ((p - /) 2 + 7 2 )((P + k) 2 + j 2 )(l 2 + ml){k 2 + ml) 
2 log 2 



" (4^)3(771,, + 7)' 

Similarly, recall ~| 0= o— >2f, and find 



LZ{V) = J dpdldk 



1 



P 4 ((P ~ 2 + 1 2 )((P + k ) 2 + 7 2 )(^ 2 + ml)(k 2 + m 2 ) 

1 r / , n 21og2 + l 1 
v{m n +7) J. 



(4vr) 3 (m 7r +7) 3 

LO{n) and LZ(n) for n > 1 are obtained similarly. Note that £(n) and LO(n) 
are finite and LZ(n) contains artificial infrared divergences, as expected. 

This exhausts our treatment of the integrals necessary to calculate Fig.l. 
A few comments are in order. In the CPS system the regularization and 
subtraction are done in the same step. The scale /i can be viewed as a label 
for the divergence and also as the renormalization scale. The counterterms 
introduced to subtract the infinities satisfy the local symmetries of the La- 
grange density. Specifically, gauge invariance and Euclidean invariance are 
easily seen to be preserved by the CPS regulator because the poles in CPS 
are isomorphic with poles in dimensional regularization, at least to the order 
we are working. This equivalence also implies that the KSW power counting 
scaling is preserved and the CPS scheme may be consistently implemented 
with other PDS calculations. Particularly, in section V the previously de- 
termined numerical values for the NLO counterterm coefficients C2 and D2 
HI are used with CPS calculated amplitudes for the deuteron quadrupole mo- 
ment. Also note that the derivation of these integrals does not change as we go 
above threshhold, so analogous results can be obtained with the prescription 
7^ — ip. The application of these methods to higher order diagrams, such 
as the 3 potential pion exchange graph analogous to Fig.l, is straightforward. 
For example, the master integral for the 3 pion exchange graph, involving 
seven propagators, is expressed in terms of Polylog of order 3 and lower. In 
general, an n-pion exchange graph can be expressed in terms of Polylog n and 
lower order polylogarithms. 



8 



.2. Photon coupled to a deuteron bound state with one pion exchange. This gives the deuteron 

quadrupole moment at NLO. 



III. FINITE MOMENTUM TRANSFER GRAPHS AND ABOVE 
THRESHHOLD SCATTERING 

Now it is shown how to calculate bound state graphs with finite momen- 
tum transfer. For pedagogical reasons we begin with a simpler graph than 
Fig.l. The analogous one pion exchange graph shown in Fig. 2 serves to illus- 
trate CPS extension to the q^O case. 

Again, in the effective field theory of 0, we find a Born amplitude of 

-3eg 2 M$j f ~ 2l i l j -l 2 S ij 



A 2/2 J dpdl (0 - q/2) 2 +7 2)(p2 + a 2 )((;p + /) 2 + a2)(/ 2 +m 2 ) - ( 30 ) 

To extract the contribution to the electric quadrupole or monopole form fac- 
tor, we multiply by the appropriate orthogonal tensor structure. These are 
Qij = HQj — and Sij, respectively. Depending on which structure is de- 
sired, different integrals will be needed. But in both cases the master integral 
is 

r ~ 1 
dpdl 



((p - q/2) 2 + 7 2 )(;p 2 + a 2 )((p + I) 2 + a 2 )(Z 2 + m 2 ) 
1 /• 3 3 exp [-(a + m^x - ay - j\x - y \ + i q ' {x ~ y ^ ] 



d 6 xd 6 y ^ _ _ — 1 1 (31) 

(47TJ J x 2 y\ x — y \ 



where the trivial delta function integrations have been performed after going 
to coordinate space. Now we write ~q-(x — y)=q\~x — y\u q and 7 = 7 — ^y 2 -, 
where u q = cos 6 q will be averaged over later. In familiar fashion, the angu- 
lar integrations are performed and a sum of 2 artificially divergent integrals 
remains : 

32^ 2 7i Wo 

Using the methods developed above, and being sure to treat the two terms in 
an equivalent manner, the following result is quickly obtained : 
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Fig.3. Box diagram contributing to the NNLO scattering amplitude. 



I = 



+ 



-)loj 



m n + a + 7 



327r 2 7 7 — a 7 + a 7 m,r + 2a 

Now we must perform the u q integration. Using eq.(23) and after some algebra 
the integral becomes 



' {29[Polylog 2 (- 



32ir 2 aq 



7 -^ 2 ) + Polylog 2 ( - a -r Z9/2 )]+lo g 



m % + 2a 



m. 



1 + 



2a 



m 



arctan 



Although this is more complicated than the zero momentum transfer integral, 
eq.(8), we still obtain a nice closed form solution. Similarly, we can derive an 
expression for the two pion exchange graph shown in Fig.l at finite momentum 
transfer. The result is in terms Polylog 3 and lower order polylogarithms and 
logarithms. 

Now we briefly comment on above threshold scattering graphs such as the 
box diagram in Fig.3, where p and q are the initial and final state momenta, 
and i and j are the spin indices for the initial and final states, respectively. The 
same techniques can be used for S^S, S^D, and D^D wave scattering. In 
any case, the momenta in the numerator of the loop integral may be decoupled 
with the propagator momenta, so that only 1, 2, and 3 propagator integrals 
are left. The master integral is 

1 



' (12 _ p2)((J _ p )2 + m 2)(( 9 _ 1)2 + m 2) • 

After the delta function integrations in coordinate space 

j. _ 1 f(f xa [ 3 y eyip ^ ipv ~ m " KX ~ m ^ ^ + y ' + i ~p ~ i ~q + 



(^) 3 J " xy\x + y\ 

Averaging over the angles of p we have 

j = 1 fd 3 xd 3 y exp ^ ipy ~ m ^ x ~ m ^\ % + y \ - 1 ~q + 



{AirfipJ ■ x 2 y\ x + y \ 

But this is exactly the same form of eq.(31), and we can proceed almost 
identically. The similarity in the derivation of these two loop integrals is no 
accident. Both for bound state graphs with finite momentum transferred into 
the loop, and for above threshhold scattering graphs, we have an external 
momentum flow into the internal loop structure. 
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IV. CALCULATING ELECTROMAGNETIC PROPERTIES IN 
EFFECTIVE NUCLEAR FIELD THEORY 

Here we only briefly review the pertinent aspects of the KSW effective 
field theory and the tools necessary to calculate electromagetic properties. 
The Lagrange density with pions and nucleons is written as 

C = C + d+ C 2 + ... 

where C n contains n-body nucleon operators. Here 

1 f 2 f 2 

C = -(E 2 - B 2 ) + i-TrD^D^ + ^-ATrm 9 (S + £t) + ... 

2 o 4 

where m q = diag(m u , m<j), m 2 = X(m u + m^), / = 132MeV is the pion decay 
constant, and the covariant derivative acting on £ is 

D^Y, = dyX- + ie[Q em , S]^. 
The one-body terms are 

A = N\iD + + !|ijvV(£D£t - ^D^iV + ... 

where the covariant derivative acting on the nucleon fields is 

D„N = [d^ + ieQ em A^)N. 
The two-body operators relevant to the present work is 

C 1 i : 2 

C 2 = -{C + D 2 XTTm g ){N T P l N)^N T P i N) + -^[{N T P t N)\N T PiD N) + h.c] 

o 

+ ^-{(^PiN^^PaiDjOi - B 2 %]iV) + h.c.} + ... 

o o 

where Pi = 4gCT2 o {T 2 projects onto spin and isospin states in the spin triplet 
channel. Also, all of the above operators are assumed to be operators in 
the spin triplet channel, which is relevant to the deuteron. Three of the 
counterterms have been determined previously |IJ to be 

C = -5.51fm 2 ,C 2 = 9.91fm 4 ,L>2 = 1.32fm 4 

at renormalization scale [i = m n . Note in eq.(43) the presence of a new local 
four nucleon operator which mixes 3 <Si and 3 D\ wave states. This does not 
appear until NNLO because it scales like 1/Q in the power counting, unlike 
the other two derivative coupling, C 2 , which scales as 1/Q 2 . The reason is 
that the C 2 operator can be renormalized by Cq bubble chains on both sides 
whereas the C 2 ^ D operator is only renormalized on the S wave side. The 
value of this counterterm will be determined in this paper. We also note that 
a direct quadrupole moment operator where the photon couples to the four 
nucleon vertex does not appear until NNNLO, so we can safely neglect it. 
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+ Higher Order Graphs 




+ Higher Order Graphs 



Fig. 4. A graphical expansion of the irreducible three-point functions (top) and the irreducible 
two-point functions (bottom). The first line in each case is the LO amplitude while the second 
lines are the NLO amplitudes. The dashed line represent potential pions, the 'X' marks are the 
deuteron interpolating fields which create or destroy a deuteron in the spin triplet channel. 

In ref.[5] the formalism necessary to calculate electromagnetic form factors 
in effective field theory were developed. The matrix element of the electro- 
magnetic current is related to the irreducible 2-point and 3-point functions 
via 

{P',j\Jem\P^) = i 

where is the irreducible 3-point function, which is defined as the sum 
of all diagrams with two nucleons interacting with an external photon and 
each other such that the graph cannot be taken apart by cutting at a Co 
vertex. The irreducible 2-point function S is defined similarly, except in this 
case there are only two nucleons interacting with each other. The derivative 
of this 2-point function in eq.(45) is the wavefunction renormalization that 
arises from introducing the deuteron interpolating field 

Vi = N T PiN, (46) 

which annihilates a deuteron in a definite spin and isospin state. Also note 
that in eq.(45) q = p' — p is the photon momenta and E and E' are the center 
of mass energies. 

The graphical expansions of the 2-point and 3-point functions up to NLO 
are shown in Fig. (4). The zeroth component of eq.(45) is related to the 
definitions of the form factors by 



r&(^,q) 



dZ(E)/dE 



R.F' 



(45) 
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(p/,j|J°JP,*) = eiFciq 2 )^ + -^Fgiq 2 )^ - -<? 2 c%)]. 

d 

For the present purposes, we need not consider the spatial components of 
this matrix element, since the quadrupole form factor Fq(q 2 ) is related to the 
quadrupole moment by 

*o(o) _ „ 



V. THE NNLO CALCULATION OF THE DEUTERON QUADRUPOLE 
MOMENT 

The deuteron quadrupole moment was calculated in ref.[5] to NLO. The 
leading order contribution vanishes and the NLO result is given by 

NLO = 9A M (&y 2 + 9 7"V + 4rr4) 
^ Q 30 7 r/ 2 (m 7r + 2 7 ) 3 

The numerical value for this is about 0.418fm 2 , which is more than %40 larger 
than the experimental value of 0.2859fm 2 . This is in very rough accordance 
with the expectation Q that each order in the perturbative expansion is sup- 
pressed by about %30. A number of recent calculations suggest that 
the expansion does in fact converge as expected in calculations up to NLO. 
However, there has not yet been a test of the theory at NNLO. Thus, it would 
be interesting to calculate the deuteron quadrupole moment at NNLO. Since 
it is only the second non-vanishing order, the errors are expected to be within 
%10, which is rivaling the accuracy of the most recent potential model calcula- 
tions H . Calculations using potential models are consistently lower than the 
experimental value of the quadrupole moment by ~%7. This suggests that 
dynamical considerations beyond potential interactions are needed to accu- 
rately describe such systems. The effective field theory approach provides 
the systematic framework to account for these effects. 

However, as mentioned previously, there is a new direct S^D wave oper- 
ator at NNLO whose value is unfixed. One may expect to extract this value 
from low energy N-N scattering data in the spin triplet channel. The diagrams 
which give rise to this AL = 2 S^D transition are shown in Fig. (5). Fitting 
these to the triplet channel mixing parameter e\ would yield the desired value 
of Cf . However, new complications arise in the NNLO above threshhold 
scattering channel that are not present at lower orders nor in NNLO bound 
state problems. Work in this area is ongoing and will be presented later. 

Alternatively, we can calculate the graphs necessary for the quadrupole 
moment and fit the new operator to the experimental value of \iq. This 
is the approach adopted in the present work. Thus, some of the questions 
raised in the preceding paragraph will have to wait until the NNLO triplet 
channel scattering analysis is complete. Several sources |l^Jll| suggest that 
perturbative potential pion exchange will give a large result compared 
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Fig.5. The top line shows the graphs contributing to the 3 Si—> 3 D\ wave scattering amplitude at 
NNLO. The bubble chain of Co operators is defined on the second line. 

with the counterterms. This expectation can be verified by the results pre- 
sented here. 

Expanding eq.(45) in powers of Q, which is the power counting scale in 
the KSW power counting scheme, we arrive at 

(p, | P , - ,\ F ?°> 1 + ; M^B-T^ {2) /dE 

(P',j|J„„|P,i)-![ (JE(i)/ii gj+4 (<ffi (1) /dB)2 J + "" (50) 

where we consider only the graphs that give a quadrupole contribution and 

1 - i (o) + i (i) + - 

£ = £ (1) + £ (2) + ... (51) 

The first term in eq.(50) is the NLO result previously calculated while the 
other two terms are NNLO corrections. A number of diagrams occur at NNLO 
in the effective field theory power counting. Here we are only concerned 
with those that give rise to a quadrupole structure Qij = qiqj — -jfq 2 , where 
q is the photon 3-momenta. These diagrams must have a virtual S^D 
and then D^S transition. The graphs contributing to are shown in 
Fig. (6). The radiation pion refers simply to the graph evaluated by taking 
the pion pole in the energy integration rather than the nucleon poles. While 
the radiation pion graph (VI) shown in Fig. (6) occurs at NNLO and has a 
quadrupole contribution, the quadrupole piece itself is at NNNLO and can 
be safely neglected here. Furthermore, graph II is found to have a vanishing 
quadrupole contribution, which can be understood as follows. Any S^D 
transition from the first pion followed by a D^S transition from the photon 
will be exactly canceled by the reverse process, namely a S^D transition from 
the photon and then a D^S transition via the second pion. In other words, 
the photon interacts symmetrically with the deuteron, whereas assymmetries 
in the electromagnetic interaction give rise to the quadrupole moment. The 
remaining four graphs give non-vanishing contributions. 

The C2 and D2 graphs are directly related to the NLO one pion exchange 
graph in Fig. (2) by 
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Fig. 6. The graphs contributing to the NNLO calculation of the deuteron quadrupole moment. The 

photon corresponds to A . 

111 ~ 2tt NLO 

r o _ mlM N D 2 (fi--f) r0 , , 

1 IV — 7T~ 1 NLO- \ oz ) 



The C2^ D graph (V) is given by 



Ty ~ 384^ • (53) 

Finally, the two potential pion exchange graph is calculated using the CPS 
regularization and renormalization scheme developed in section II and deter- 
mined to be 



_ 9e g 4 M 4 



12807 9 + 35687 s "V + 4272 7 7 m 7r 2 + 32967 6 m 7r 3 + 1664 7 5 m 7r 4 + 3497 4 m 7r 5 + 33j 3 m w 6 



+ 607 2 m 7r 7 + 187m/ - 2/u(5207 8 + 17007 ? m 7r + 23067 6 m 7r 2 + 20077 5 m 7r 3 + 16377 4 m 7r 4 + 12157 3 m 7 

1 u?{i' i 4- 3m 2 ) 

+ 6037 2 m 7r 6 + 1627m 7r 7 + 18m 7r 8 ) /(294912vr 3 7 5 (7 + m 7r ) 2 (27 + m 7r ) 3 ) + P ' 71 



5 



/^ lo S (77 2 + ^ 2 )log 2 ^ m, 2 (3m, 2 - 4 7 2 )(Polylog 2 



24576vr 3 7 5 

^ 12 > 



27+m-n- 



3072vr 3 7 2 24576vr 3 7 3 6144vr 3 7 5 

(I8O7 5 + 3407 4 m 7r + 4097 3 m 7r 2 + 12 7 2 m 7r 3 - 2137m/ - 72m/) log 2 ^+^ } 1 
+ 737287r 3 7 4 (2 7 + m^) 2 J (M) 

The expressions for the two-point functions are reproduced here and were 
calculated in H to be 
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B 87T7 



E=-B % 



The partial contributions to the quadrupole moment by graphs I, III, IV, and 
V are 

/j,q = -lG^m 2 ,^ 11 = 0.079fm 2 ,^ y = -0.096fm 2 , 

and 
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Also, the last term of eq.(50) gives an additional wavefunction renormalization 
contribution of 

_ dZ (2) /dE 
= 0.203fm 2 . 

The renormalization scale /j, has been set to m n but it should be noted that 
the running of the counterterms cancels the scale dependence to the order 
we are working. Fixing the counterterm Cf^ D to reproduce the physical 
quadrupole moment yields a value of 

C^ D = 5.1fm 4 . 



CONCLUSIONS 

We have shown that otherwise comlicated loop diagrams arising in non- 
relativistic field theories, particularly the EFT of Kaplan, Savage, and Wise, 
can be easily calculated in a new regularization and renormalization scheme 
called CPS. The divergences are dealt with consistently by transforming 
them into divergences in a parameter with units of energy. This is done af- 
ter the angular integrations in configuration space by inserting the identity 
operator in the form J^da-^, where a appears in the exponential. For artifi- 
cial divergences such as the artificial infrared divergences arising in Fig.l, the 
only requirement is consistency in dealing with the infinities, since they cancel 
in the end. For real divergences, the CPS method precisely reproduces the 
results of dimensional regularization and PDS (or MS) up to NLO. Higher 
order loop integrals evaluated in CPS also map onto poles in dimensional reg- 
ularization. This gives one the power of the PDS subtraction scheme and the 
associated KSW power counting, without the calculational difficulties associ- 
ated with dimensional regularization and Feynman parameters. As is the case 
with PDS, we could ignore the power law divergences and keep only the four 
dimensional logarithmic divergences, yielding a subtraction scheme similar to 
MS. To illustrate the methods, the quadrupole moment of the deuteron is 
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calculated to three loops, thus fixing the numerical value of the direct four 
nucleon 3 Si^ 3 -Di operator that appears at this order. 

Although the emphasis in this paper has been effective field theories in 
nuclear physics, the techniques developed are applicable in other field theories, 
such as non-relativistic quantum chromodynamics (NRQCD). 
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